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ON  CERTAIN  UNSTEADY  MOTIONS  OF  A  CCMPRESSIBLE  FLUID 
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kj#  exact  solutions  of  the  equations 

afiJS  O  A-W. 

for  one-dimensional  unsteady  motion  of  a  compressible  fluid/in  the  case  of 


? — ’ 

plane  waves  and  motions  with  cylindrical  and  spherical  symmetry.  The 
solutions  are  determined  by  meane-^T  simple  procedures  based  on  considera¬ 


tions  of  the  theory  of  dimensions | 


considered-as - 


a  general  method  that  permits  us  to  find  in  many  questions  of  physics 
families -  of  spltrbiona  vrfiich  “ depend  on  arbitrary -parameters'.'^  In  concrete 


examples  it  is  usualiy  easy  to  Indicate  the  arrangements  of  the  problems 

and  the  hypotheses  leading  to  solutions  of  the  type  obtained.  We  shall 

j  f  '-Jld-- 

apply  .this  general  mettiod^to  both  linear  and  nonlinear  equations  with 

a 

partial  derivatives  in  the  most  varied  problems  of  physics  and  mechanics.^ 
The  method  of  Fovirier, which  is  applicable  to  certain  linear  problems,  ■ 
is  widely  known.  Sometimes  it  is  possible  to  utilize  partisJ.  solutions, 
which  are  obtained  by  means  of  the  general  method  under  consideration,  as 
well  as  partial  solutions  obtained  by  the  method  of  Fourier,  for  setting 
up  the  general  solution  of  equations  of  motion  emd  for  the  solution  of 
boundary  problems. 


^A  short  report  on  these  solutions  was  published  in  a  note  by  the 
author  [ij  . 
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.iXl  Let  U8  consider  the  one-dimensional  unsteady  motions  of  a 
compressed  fluid.  Let  us  take  time  t  and  the  linear  coordinate  r  as 
Independent  variables.  As  basic  unknown  quantities  let  us  take  the 
velocity  v  of  the  particles  of  the  fluid,  the  density  p  and  the  pressure 
p.  The  dimensions  of  density  and  pressure  depend  on  the  symbol  of  unit 
of  mass;  therefore  the  functions  p(r,t}  and  p(r,t)  must  contain  certain 
dimensional  physical  constants.  These  constants  may  be  a  part  of  the 
equations  of  motion,  which  the  vinknown  functions  satisfy,  or  the  additional 
conditions  by  which  the  solution  is  determined  (characteristics  of  the 
field  occupied  by  the  moving  fluid,  boundary  conditions,  initial  conditi(»i8, 
etc. ) . 

Let  us  consider  the  cases  where  the  unknown  functions  contain  dimensional 
constants,  among  which  there  are  only  one  or  two  constants  with  independent 
dimensions.  Special  motions  of  such  a  kind  may  represent  independent  in¬ 
terest  or  may  prove  to  be  finite  motions,  if,  in  the  limit,  among  the 
assigned  physical  characteristics  only  one  or  only  two  characteristics  with 
independent  dimensions  are  essential. 

Let  us  assume  in  the  beginning  that  there  is  only  one  dimensional  con¬ 
stant  which  we  denote  by  a.  It  is  obviowS  that  the  dimensions  of  the 
constant  a  must  contain  the  symbol  of  mass;  therefore,  without  loss  of 
generality  we  can  consider  that  the  dimensional  formula  for  a  is  of  the 
form 

k  B 
a  -  M  L^  T  , 

where  k  and  »  are  constants.  For  the  quantities  v,  p ,  and  p  wo  can  write 


the  formulas 
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V  -  r  V  ,  P  - 


_eR- 


aP 


(1.1) 


rk.lt.*2 

It  is  obvious  that  the  abstract  quantities  V,R,  and  P  canot  depend 

O 

on  r  and  t  and,  consequently,  can  only  be  constants.  Consequently,  the 

hypothesis  that  only  one  constant  a  is  essential  makes  it  possible,  with¬ 
out  resorting  to  the  equations  of  motion,  to  completely  determine  the 
dependence  of  v,  p ,  and  p  on  r  and  t.  If  the  equations  of  motion  do  not 
contain  physical  constants  or  contain  constants  with  dimensions  expressible 


in  terns  of  the  dimension  a,  then  these  equations  must  have  solutions  of 
the  kind  (l.l)  when  V,  R,  and  P  are  constants.  It  is  not  difficult  to  see 
that  in  this  case  whatever  the  constants  k  and  s  may  be,  the  motion  determin¬ 
ed  by  the  formulas  (l.l)  does  not  contain  the  condition  for  vdiich  density 
and  pressure  are  different  from  zero  emd  at  the  same  time  constant  through¬ 
out  the  entire  mass  of  fluid. 

If,  besides  the  constant  a,  there  is  still  a  second  constant  ^  then, 
without  loss  of  generality,  we  can  always  assume  that  th^  dimensions  of  b 
do  not  contain  the  symbol  of  mass  and,  consequently,  a  formula  of  the  form 

|bj  -L®T", 


where  m  and  n  are  constants,  is  true;  it  is  obvious  that  %dien  n^O  only  the 
ratio  m/n  is  essential. 

From  the  determining  quantities  jr,  t,  ^  and  b  we  can  form  only  one 
independent  abstract  variable  combination 

X  -  b  r"®t~". 


^Concerning  general  methods  of  the  theory  of  dimensions  see,  for  example, 
the  author*  s  book  . 
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Turning  to  formula  (1*1)  we  conclude  that  the  abatract  quantities 
V,  R,  and  P  can  be  functions  of  a  single  variable  parameter  Further, 
we  shall  show  that  we  can  construct  solutions  of  the  type  under  consideration 
in  which  the  fvmctions  V,  R,  and  P  or  their  derivatives  are  discontijmous. 

The  hypothesis  that  only  two  constaints  a  and  b  with  independent 
dimensions  are  essential  permits  us  to  determine  the  general  type  of 
distributions  of  the  values  v^,  p^,  and  p^,  when  t  ••  0.  By  means  of  the 
theory  of  dimensions  |2j  we  find 


-  i  s  +  1  - 1  m  -k-3 

v^  ■  b  r”  »  r  "  ,  p 


fc  _  2  ms  ^  ^  2m 

n  "  Dr  n  "  "  ^  gj 


Vq  -  b  r  >  Pq  ■  ^  ^ 

where  0(2,  Q!^  are  abstract  quantities  that  we  may  consider  as  functions 
of  r,  having  constant  values  when  r^O  and  r/coand  having  singular  points 
when  r  -  0  and  r  ■  oo  in  the  general  case.  The  quantities  a a2»  anda^# 
in  particular,  may  be  equal  to  zero  or  infinity. 

If,  when  t  ••  0,  pressure  and  density  are  constant,  finite,  and  different 
from  zero,  then  the  following  relations  must  hold: 


-  -  -  1,  k  +  8 


•  -  5,  ( 


(1.3) 


from  which  it  follows  that  the  magnitude  of  the  Initial  velocity  is  the  same 
for  all  particles  of  the  fluids.  Discontinuity  of  velocity  can  take  place 
only  when  r  ®  0  or  when  r  =  <*>  . 

Let  us  note  certain  problems  for  which  the  solution  can  be  constructed 
by  means  of  solutions  of  type  (l.l)  when 


k  “  -  3, 


s  -  0, 


5  -  -1. 


(1.4) 
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A.  Classical  problem  of  ari  **  explosion "  along  a  plane.  It  is  necessary 
to  determine  the  adiabatic  motion  of  a  gas,  if  at  the  initial  moment  of  time 
the  velocity,  density,  «ind  pressure  have  a  discontinuity  along  the  pleuie 
X  r  =  0,  while  on  the  rigjit  hand  side  for  r  >  0  we  have  the  assigned 
constant  values  v^,  0^,  and  p^,  and  on  the  left  hand  side  for  r  <  0 
we  have  the  assigned  constant  values  Vg,  Pg,  and  Pg. 

As  is  known,  equations  of  motion  of  a  perfectly  ideal  gas  do  not  con¬ 
tain  dimensional  constants;  thei'efore,  in  this  case  the  motion  depends  only 
on  the  six  dimensional  constants  v^,  p^,  p^,  Vg,  pg,  and  Pg. 

Among  these  constants  there  are  only  two  with  independent  dimensions. 

We  can  take  a  and  b  as  constamta : 

a  =*  p^  ^ =  ML  ,  k  —  -5,  B  ^ 

(1.5) 

»  "  !!i  ([»>]-  m  .  2,  n  .  -2,  a  -  -x). 

Pj 

Relations  (l.^)  are  satisfied.  The  solution  of  the  assigned  problem 
can  be  obtained  by  means  of  a  combination  of  solutions  of  type  (l.l).  Con¬ 
sequently,  in  this  case  we  can  say  in  advance  that  the  foraulas  which  give 

expressions  for  v,  p,  and  p  in  terms  of  r  and  t  can  depend  only  on  the 
2/2 

parameter  bt  /r  *  k.  Ibrmulas  expressing  v,  p,  and  p  in  terms  of  r  and 
t  can  change  their  form  when  passing  throu^  discontinuities^  (weak, 
stationary,  emd  strong  Jun^js). 

^For  an  examination  of  possible  cases  of  motion  in  this  problem,  see  the 
work  of  N.  E.  Kochin  [^3j- 
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B.  Problem  of  the  aotion  of  a  gas  In  a  long  cyllndrlo*!  tube,  cloaed  »t 

one  end  with  a  pieton.  At  the  first  moment  of  time  the  gas  is  at  rest,  and 
the  piston  begins  to  moTO  suddenly  with  constant  velocity  U.  By  reversing 
this  motion,  wo  arrive  at  the  problem  of  the  sudden  introduction  of  a  bander 
in  front  of  an  advancing  homogeneous  stream  of  gas.  ’ 

In  this  problem  the  motion  is  determined  by  the  following  dimensional 
constants:  the  Initial  density  the  initial  pressure  p^,  and  the  velocity 
of  the  piston  U.  We  again  obtain  the  case  vAiere  the  solution  depends  only- 
on  two  constants  with  independent  dimensions.  As  is  known,  if  U  is  directed 
toward  the  gas,  then  the  solution  of  the  problem  has  the  simplest  form  -  - 
in  front  of  the  piston  there  is  formed  a  sone  with  constant  pressure,  density,  , 
and  velocity  of  gas  particles  equal  to  U.  This  zone  is  separated  from  the 
quiescent  gas  by  a  strong  discontinuity,  expanding  with  constant  velocity. 

C.  Problem  of  the  detonation  of  a  gas  in  a  cylindrical  tube  (detonation 
begins  either  in  a  plane  cross  section  of  the  tube  and  propagates  in  both 
directions  in  a  quiescent  medium,  or  the  detonation  begins  from  the  end  of 
the  tube  closed  by  an  immovable  piston) . 

During  detonation  there  is  added  a  constant  that  defines  the  resulting 
specific  heat  of  chemical  reaction  q^^  -  ^h*  constant  q^  -  q^  in  mechanical 

units  of  'measurement  has  the  dimensions  of  the  square  of  velocity,  vrtiich 
depend  on  the  dimensions  of  pressure  and  density.  Consequently,  in  this 
case  we  have  only  two  dimensions^,  constants  with  independent  dimensions, 
and  the  solution  can  be  formed  from  solutions  of  the  typo  (1.1). 
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The  solution  can  be  generalized  by  assuming  that  the  piston  at  the 
initial  moment  of  time  begins  to  move  with  a  constant  velocity  U,  or  by 
assuming  that  Instead  of  a  movable  piston  near  the  end  of  the  tube,  there 
is  created  a  constant  pressure  that  is  maintained  with  time  and  differs 
from  the  pressure  in  the  quiescent  gas. 

For  plane  waves  the  solution  of  the  problems  enxjmerated  is  elementary. 
We  shall  examine  the  generalization  of  the  specified  problems  in  the  case 
of  cylindrical  and  spherical  symmetry. 

For  the  sake  of  concreteness  we  shall  further  restrict  ouiselves  to 
a  detailed  study  of  the  corresponding  problems  with  spherical  symmetiy. 
Before  that,  we  shall  consider  the  following  problemi  At  the  moment  ’ 
t  =  0  inside  a  sphere  3^  of  radius  r^,  the  gas  is  quiet  and  has  a  given 
density  P2  and  pressure  P2«  outside  the  sphere  3^  there  is  the  same  gas, 
in  which  density,  pressure,  and  velocity  are  constant  and  equal  to 
Pi,  and  Vj^.  It  is  required  to  determine  the  starting  motion  of  the 

It  is  obvious  that  in  this  case  the  solution  depends  on  six  dimensicmal 
constants  r^,  Pg,  p^,  Pj^,  p^*  which  there  are  three  constants 

with  independent  dimensions.  Since  the  radius  r^  in  the  general  case  has 
an  essential  value,  the  solution  of  this  problem  does  not  have  the  form 
(1.1).  If  we  proceed  to  the  limit  as  r^->0,  then  as  a  result  we  shall 
obtain  a  problem  with  two  independent  dimensional  constants,  iirtiich  can  be 
solved  by  means  of  solutions  of  the  type  (l.l). 

Thus  we  arrive  at  the  following  problems  with  spherical  symmetry  for 
iriiich  we  shall  give  the  methods  of  solution. 
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1«  All  particles  of  homogeneous  gas  have  th«  game  velocity  directed 
toward  tho  center  0.  <  0  (foeussed  at  >  point)  or  from  the  center 

v^<  Q  (dispersion  frc«n  center) . 

In  both  cases  near  the  point  0  there  is  formed  a  spherical  expanding 
region  of  quiescent  fluid,  in  the  first  case  with  increased  prcssxire,  in 
the  second  with  reduced  pressure.  Frca  the  following  it  will  be  clear 
that  there  exists  a  critical  velocity  v*  >  0.  If  >  v*,  then  a  vacuum 
forms  at  the  center  0. 

2°  Generalization  of  problem  B  concerning  the  piston.  When  t  ”  0 

the  gas  is  quiescent,  but  the  particles  of  the  sphere  of  infinitesimal 
radius  with  center  at  the  point  0  suddenly  obtain  a  constant  velocity  U 
directed  radially  from  the  point  0.  In  the  following,  when  t^O  these 
particles  continue  to  move  with  the  same  constant  velocity  U.  The  fl\iid, 
as  it  were,  is  laoved  apart  by  the  sphere  S^,  the  radius  of  which  increases 
from  zero  proportionally  to  time.  By  means  of  the  theory  of  dimensions  it 
is  not  difficult  to  see  that  on  the  surface  of  the  sphere  the  pressure 

f 

must  be  constant. 

The  corresponding  pressure  is  determined  from  the  function  of  the 
velocity  D  and  of  the  characteristics  of  the  quiescent  gas.  Consequently, 
this  problem  coincides  with  the  problem  of  the  expansion  of  the  quiescent 
mass  of  fluid  by  the  assigned  constant  pres8\ire  >diich  arises  at  the  point  0. 

3?  The  problem  of  spherical  detonation.  When  t  ■  0  at  a  certain  point 
0  in  the  gas,  there  arises  a  detonation  that  then  propagates  in  the  quiescent 
gas.  In  this  case,  az*ound  the  point  0  an  expanding  sphere  forms  from 
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the  quiescent  gas.  The  density  and  pressure  inside  S.,  are  constant  and 
differ  from  the  density  and  pressure  of  the  quiescent  gas  through  which 
the  detonation  wave^  has  not  yet  passed. 

2.  To  obtain  the  conditions  that  must  satisfy  the  values  which  enter 
formulaB  (1.1),  it  is  necessary  to  resort  to  the  equations  of  motion.  It 
is  obvious  that  to  Justify  the  accepted  hypotheses,  it  is  necessary  that 
the  equation  pf  motion  should  not  contain  essential  constants  with  dimensions 
independent  of  the  dimensions  a  and  b. 

For  the  saJke  of  definiteness  of  the  problem,  let  us  take  the  equations 
of  Biotion  of  a  fl\iid  in  the  following  form: 

^  +  V  =  0,  (2.2) 

where  7  is  some  abstract  constant.  When  V  ■  1  we  have  plane  waves)  when 
V-  2,  cylindrical  waves;  and  when  V  "  3»  spherical  waves. 

When  7  *  ^p/^Y*  equations  (2.1)  and  (2.2)  can  be  considered  as  the 
equations  of  the  adiabatic  motions  of  a  gM,  where  entropy  can  be  different 
for  different  particles  of  the  gas.  Equations  (2.1)  and  (2.2)  do  not  contain 
dimensional  constants;  therefore  these  equations  must  have  a  solution  of 
the  type  (l.l),  where  there  is  only  one  constant  ^  and  also  in  the  oase 
where  there  are  two  essential  constants  a  and  b. 

^he  problem  of  spherical  detonation  was  solved  by  Zeldovichp^.  Like 
Problen  C  on  detonations  in  a  cylindrical  tube,  the  problem  ox  spherical 
detonation  can  be  complicated  by  the  addition  of  the  contitions  of  Problem 
2°. 
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Substituting  forniulaB  (l.l)  into  eqviations  (2.l)  and  (2.2),  ve 
obtain 


(a  +  mV)V»  +  “  I J  =  V®  -  V  -  (k+l)  |  , 

\  [mV*  +  (n+mV)  |'J  =  -  a  -(k-V+3)V,  (2.3) 


l(n+mV) 


P* 


R« 


P-"  R 


=  -  8(1-7)  -2  -  k(l-7)  +1-37  V. 


If  the  constant  ^  it  absent,  then  the  left-hand  i&cinbers  of  equations 
(2.3)  niust  be  replaced  by  zeros,  as  a  resxalt  of  vhich  three  slntple  finite 
relations  are  obtained  ecmnectlng  V,  P,  R,  k,  and  s: 


^  =  g4i^7-ir-  ** 


(k+l) 


2(v-k-3 


k-3) 

FI) 


(2.M 


In  this  case  formulas  (l.l)  give  a  fauaily  of  exact  solutions  of 
equations  (2.1)  and  (2.2),  depending  on  the  tvo  arbitrary  constants  k 
and  Replacing  t  by  t  -  t^  wa  obtain  a  solution  that  depends  on  three 


constants. 

When  y  =  1  In  formula  (l.l),  the  coordinate  r  can  be  replaced  by 
r  -  r^,  after  which  we  obtain  a  solution  that  depends  on  four  arbitrary 
constants. 

For  an  ideal  gas,  the  entropy  is  expressed  by  the  combination 


1.7  P  j^.k(7-l)  -*>37-1  ^8 (7-1) -2 

a  g7  7 


which  in  the  general  case  is  different  for  different  particles  of  fluid. 


If 


k  = 


37-1 


7-1 


f 


s  = 


2 

7-1  ' 


then  for  the  insulting  family  of  solutions  that  depend  on  one  constant, 
the  entropy  is  identical  for  all  particles  of  the  "fluid. 
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When  mjfcOf  integration  of  the  system  of  eqmtions  (2.5)  can  toe 
reduced  to  an  integration  of  a  single  ordinary  differential  equation 
of  the  first  order. 


dz 


z  {  [2(V-l)+T<r-l)v](V-q)^-(r-l)V(V-l)(V-q)-[2(V-l)+x(7-l)j  zj 

(V-q)  [v(V-l)(V-q)+(x-W)zJ 


5) 


where  the  following  designations  are  used: 


z  » 


=  _  2.  •  -  s-i-2»q(k+l) 


After  equation  (2.5)  haa  toeen  integrated,  the  dependence  of  V,  R 
on  \  will  toe  determined  toy  means  of  quadratures  from  the  relations 


d  In  (V-q)^-z 

dV  "  V(V-l)  (V-qi+(x-i;V)z" 


(2.6) 


fv_q)  -  rs+(k-v+3)vj .  (2.7) 

d  In  z-(v-q)^ 

In  the  general  case  the  integration  of  equation  (2.5)  must  toe  perform¬ 
ed  approximately.  When  m  =■  0,  i.e.,  when  q  *  oO ,  eqviation  (2.5)  is  easily 
integrated,  and  we  obtain  the  Integral 


2^  =  z  7AV^(V-1)^^’''^^  (2.8) 

where  A  is  a  constant  of  integration.  Assuming  n  =  1  (i.e.,  \  =  to/t, 
where  to  is  the  characteristic  time),  we  find  from  equation  (2.3)  that 


d  In  \  = 


dV 


V[y-l-A(k+l)7(V-l?^^^"^^] 


(2.9) 


d  In  R  \ 


8 


(k-»3-z^)dV 

V-1-A(k+1)  V(V-l)*^^^"^^ 


(2.10) 
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In  this  case  the  solution  of  the  equations  of  motion  of  the  fluid 
(2.1)  and  (2.2)  can  be  written  in  the  form 


"t  r 


V  [(V-l)-A(k+l)V(V-l)^ 


P  =•  ^77  CJqp  (k+5-i^)  f - - yr  va},  (2.1 

^  y  (V-l)-A(k+l)V(V-lf 

p  =  ^  AV^(V-lf^’"^V,  V  -  ^  V  . 

These  formulas  give  an  exact  solution  of  equations  (2.1)  and  (2.2) 
which  depend  on  the  four  eubitrary  independent  constants  ^  sB  and  k. 

For  80B»  particxilar  values  of  7  the  integrals  entering  these  fora\xlas 
can  be  calcu].ated  in  a  simple  form. 

If  k  =■  -1,  then  for  any  7  the  integrals  are  easily  calculated,  after 
which  we  obtain  the  following  exact  solution  of  the  equations  (2.1)  and 


(2.11) 


(2.12) 


(2.2) 


t-t  " 
o 


(t-t  ) 
*  '  o 
P  - 0“ 


2-  7/ 


p  »  P 


(t-t  ' 


(2.15) 


where  t^,o<^  and  p  are  arbitrary  constants.  In  the  motions  described  by 
the  formula  (2.15)  pressure  depends  on  time,  but  is  constant  throughout 
the  entire  mass  of  fluid.  All  particles  move  viniformly  and  rectilinearly, 
but  with  different  velocities.  For  each  particle  of  the  fluid  the  entropy 
is  constant,  but  is  different  for  different  particles.  Where  the  symmetry 
is  cylindrical  the  distribution  of  density  is  stationary.  A  more  general 
solution  is  obtained  if  in  the  solution  (2.15)  the  constant oj is  replaced 
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Solution  (2.13}  eolncldea  In  form  with  the  aolutlon  fumlBhed  by 
formulas  ^l.l)  vhen  ^  R  and  P  are  constants,  determined  by  formulas 
(2.4),  but  differs  from  it  by  the  fact  that  it  depends  on  the  two 
dimensional  constants  o<.snd  ^  \diich  In  the  given  case  enter  in  a  special 
form. 

When  V  =  1  in  formulas  (2.13),  the  coordinate  r  can  be  replaced  ty 
r-r^,  after  which  we  obtain  a  solution  depending  on  four  arbitrary  con- 
stamts)  this  solution  differs  from  the  corresponding  solution  that 
satisfies  formula  (2.4). 

Solutions  in  the  form  (2.l)  can  be  sou^t  when  there  are  surfaces 
of  weak  or  sti^mag  discontinuities.  These  characteristic  surfaces  will 
satisfy  certain  values  of  the  coordinate  r  »  r*,  which  depends  on  time. 

The  dynamic  conditions  connecting  the  values  v^,  p^  and 
Vg#  02#  Pg  °°  different  side  a  of  surfswses  of  strong  Junaps  have  the  forh: 

Condition  of  constancy  of  mass 

Pl(v^-c)  =  Pg(v2-c),  (3»l) 

where  c  is  the  velocity  of  displacement  of  the  iunpj  and 

Condition  of  conservation  of  mosaentum 


irtiere  F  is  an  aihitrary  function. 


P2(v2-c)  "  ^2"^!  • 

We  shall  write  the  equation  of  energy  assuming  that  the  gas  is 
ideal: 


’'1  Pi 


*  I  '3-3> 


ii  -  ■’2  *  ■  (rj  -  iypV 
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whare  •  qg  denotcB  the  energy  referred  to  a  unit  of  mago,  t^ich  can 
be  generated  by  the  chemical  reeu:tions  acconpanylng  the  Jung).  It  is 
necessary  to  introduce  the  qiiantity  by  considering  the  phenomenon 

of  detonation.  We  shall  assune  that  the  difference  Is  expressed 

In  mechanical  units  of  measurement  and  therefore  the  dimensions  of 
qi  -  qg  are  equal  to  the  square  of  the  velocity.  The  constant  7  may  be 
different  on  different  sides  of  the  stnrface  of  discontinuity. 

In  the  Junps  t,  a  and  b  vlU  be  the  dimensional  determining  parameters 
from  which  It  ls<  ijs^sslble  to  form  a  dimensionless  combination;  there¬ 
fore,  for  solutions  of  (l.l)  on  the  Junps,  the  relations 

_  i  i  _  E 

X  =  :=  const,  r*  =  ^  ^  (3»^) 

will  hold. 


Since  c  »  dr*/dt,  then 


c  = 


--  £ 
m  t 


(3.5) 


We  can  further  assume 

q^  “  f~;  /  (3.6) 

where  X  and  q  are  abstract  constants. 

o  * 

Solutions  of  the  fcnm  considered  can  exist  during  the  generation  of 


heat  by  the  chemical  reactions.  If  this  generation  can  be  subordinated  to 
a  law  satisfying  the  formula  (3.6). 

■  Motion  where  there  are  surfaces  of  discontinuity  is  Impossible,  if 
there  is  only  one  essential  dimensional  constant  a. 

As  pointed  out  above  In  (1.3),  if  motion  starts  from  a  state  of  rest 
In  which  the  fluid  is  homogeneous,  then  m/n  =  -1.  Consequently,  for  motions 
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of  the  type  considered  that  starts  from  the  state  of  rest,  in  vhich 
and  are  not  equal  to  zero  or  infinity,  it  follows  that  the  velocity 
of  the  shock  waves  is  constant  and  the  specific  heat  of  the  reaction 
must  also  be  a  constemt. 

Substituting  forinilas  (l.l)  in  the  relations  (3-l),  (3*2),  and  (3*3) 
and  Introducing  the  designations 


v-c 

nr 

r  /t 


ve  obtain 


<l  + 


> 

(3.8) 

*1 

’'l^ 

(3.9) 

2 

2 

^2 

'^2 

(3.10) 

2 

7«-2 

2 

From  equations  (3*9)  aad  (3*10)  it  is  easy  to  esipress  Zg  and  u^  in 
terms  of  z-^,  u^  and  q. 

For  shock  waves  without  chemical  reactions  "  7g  =  7  and  q  =  O), 
two  solutions  of  equations  {3«9)  and  (3.10)  have  the  simple  form 


'"l  “  “2^  *2  = 


U2  =  u^ 


1  + 


*1"  ^ 


r+1 


r_  2 

1 - 

i 

(3.11) 


(5.12) 
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Since  equations  (3.9)  and  (3.10)  are  syinmetrleal  with  respect  to 
z^,  in  terms  of  Zg,  Vg  it  is  stiTflclent  to  transpose  the  Indices  1 
and  2  in  formulas  (3.11)  and  (3.12)^. 

In  the  following,  we  shall  assuine  that  the  index  1  corresponds  to 
that  side  of  the  to  ^Ich  we  cross  during  the  growth  of  r. 

For  a  general  analysis,  another  relation  is  useful: 


-  Ug^  -  -(z3_-U3_2)  /  1  + 


¥} 


(3.13) 


The  value  of  z  in  its  physical  sense  is  essentially  positive. 
According  to  the  theorem  of  Tsenqslen,  the  inequality  "  ^2  ^  must 

be  satisfied,  \dilch  yields 

Ui  -  -  V,  -V,  <  0.  (3.11.) 

On  the  half -surface  -  aO<V<+fl(),r>0  (Fig.  1)  are  marked  the 


Fig.  I 


q 

^In  the  preceding  consideatlons,  formulas  (2.1)  are  utilized  for  formulat- 
ii:g  conditions  in  the  Jvuzq)  in  the  dimensionless  form  (k.8),  (4.9),  and 
(4.10)  on  page  25j  therefore  the  relations  mentioned  are  true  for  any 
motions  of  a  gas  with  straight  JmqpB.  For  8i>ecial  motions  of  the  type 
considered,  we  have  u  «  V  +  n/m. 
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regions  that  transform  into  each  other  by  means  of  relations  (3*11) 
and  (3.12). 

The  straight  line  =>  0  transforms  into  the  parabola  z  =  27(V+n/n^/(y-l) . 
The  points  of  the  parabola  z  =  (V+n/m)®  will  transform  into  theoselves. 

From  the  theorem  of  Tsein>len  it  follovs  that  the  points  s^, 
must  be  located  in  the  regions  that  are  hatched  vertically^  the 
corresponding  values  of  Zg  and  7^  will  be  found  in  the  regions  that 
are  hatched  horizontally.  Possible  transitions  from  z^,  to  z^,  Vg 
are  Indicated  by  corovs. 

Points  of  the  type  M,  located  above  the  i>erabola  z  •  27(V+n/m)^/  (7-I), 
cannot  correspond  to  the  boundaries  of  the  J\ai^>s. 

Let  us  now  consider  in  more  detail  another  case  vhere  the  condition 
of  rest  is  a  particular  solutlcm  in  the  family  under  coxislderation.  In 
this  case 


whezre  a  is  the  velocity  of  sound. 

If  z  >  (V  -  1)^,  then  the  velocity  of  the  fluid  particles  with 

respect  to  the  Jump  is  less  than  the  velocity  of  sound;  when  z  <  (V  -l)®i 

2 

it  is  greater  than  the  velocity  of  sound.  On  the  paorabola  z  =  (V  -  l) 
the  relative  velocity  of  the  particles  is  equal  to  the  velocity  of 
sound. 

If  V  >=  0  when  r  jj^  0  and  than  V  =  0;  whence  it  follows  that 

the  point 0  on  the  z-euc^  correspond  to  the  conditions  of  rest. 
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On  the  straight  11ns  V  =  1  ve  have  v  =  r/t.  This  means  that  the 
velcxilty  of  the  fluid  particles  Is  sqruxl  to  the  velocity  of  transition 
of  the  conditions  of  motion  X  =  constant  that  satisfy  the  points  of 
this  sbraight  line.  Consequently,  if  there  are  stationary  discontinuities 
in  the  stream,  then  they  must  satisfy  the  points  of  the  straight  line 
V  =  1. 

Let  us  assxane  that  ve  have  a  condition  of  rest  on  side  1}  this 

condition  can  change  into  motion  by  a  Jiunp  only  for  points  of  the 

2  2 

z-axis  located  on  the  segment  O^A,  for  •which  =  a^^  /c  <1.  The 
points  of  the  Segment  transform  by  the  jun^)  into  the  points  of  the 
arc  of  the  parabola  AOg,  the  eq\iation  of  which  has  the  form 

Zg  ”  (l-Vg)  (1+  Vg)  .  (3.15) 

The  point  Og  for  which  Xg  =  27(7-1)/ (r  +  l)^,  Vg  =  2/(7  l) 
corresponds  to  the  point  0^^®!  =  ^3^  =  0).  If  the  velocity  of  sound 

in  a  quiescent  medium  is  finite,  then  the  limiting  case  when  the  speed 
of  propagation  of  the  Junq)  is  equal  to  infinity  corresponds  to  the 
Jxm?)  from  0^^  to  the  point  Og.  The  velocity  of  propagation  that  is 
equal  to  the  -velocity  of  sound  corresponds  to  the  point  A  (z  =  1,  V  =  O). 
When  the  point  z^^,  approaches  the  point  A,  the  Jumps  in  velocity, 

density,  and  pressure  converge  toward  zero.  Only  a  wsak  discontinuity 
may  correspond  to  the  point  A.  By  the  position  of  the  point  A  on  the 
segment  O^A  the  velocity  of  propagation  of  the  Jxmp  into  the  quiescent 


fluid  is  determined. 


In  the  case  of  the  propagation  of  a  direct  condensation  Juiqp  in 
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a  quiescent  gas  =  0  (vhen  ^  ® 


^2 


(3.1),  (3.2),  and  (3.3)  can  be  written  in  the  form 


7),  the  relations 

(3.16) 

(3. IT) 

(3.18) 


These  formulas  are  convenient  for  evaluating  the  elements  of  a 


Jump  by  means  of  the  velocity  of  propagation  of  the  Jump  or  by  means  of 


the  velocity  of  the  particles  of  the  gas  behind  the  Jump  Vg.  For  a  better 
Illustration,  the  relations  represented  by  tlie  formulas  (3.16),  (3.17)# 
and  (3.18)  are  shown  graphically  by  solid  lines  in  Fig.  2. 


Fig.  2 


T-57 

2-5-56 

-  20  - 


If  the  ratio  ^  aaaigned,  then  In  order  to  determine  the 

ratio  a^/c  it  ie  sufficient  to  draw  a  atralght  line  with  a  slope  equal 
to  ■v’g/a^  throu^  the  origin  of  coordinates  on  the  lower  graph.  The 
Intersection  of  this  etralght  line  with  the  parabola  (3*16)  determines 
the  ratio  Sj^/c. 

We  shall  now  examine  another  ease,  when  q^^  -  q^  ^0  and  7^  7^,  i 

and  the  shock  wave  propagates  In  a  quiescent  medium  =»  0  or  =  -1. 

In  this  case  equations  (3*9)  and  (3>10)  asstune  the  form 


(3.19) 

(3.20) 


2  2 

When  >  9^  “  ^2^  constants  in  the  z,  u  plane,  these 

equations  are  represented  by  parabolas  located  as  shown  in  Fig.  3,  if 


(3.21) 


The  equality  sign  corresponds  to  the  case  where  the  parabolas  Inter¬ 
sect  at  the  point  D.  If  the  inequality  (3.2l)  is  not  satisfied,  then  the 
system  (3.19)  end  (3.20)  has  two  solutions  -  for  one  of  them  Vg  <  0,  and 
for  the  other  1  >  Vg  >  O. 

For  the  existence  of  two  real  dissimilar  solutions  of  the  equations 
(3.19)  and  (3.20),  tlie  constants  a^,  c,  and  0.3^  “  <l2  ““sli  satisfy  the 
relation 
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> 


The  roots  of  equations  (3*19)  and  (3»20)  converge  when  in  the 
relation  (3 -22)  the  sign  of  the  Ineqiiality  is  changed  by  the  sign  of 
eqxiality}  in  this  case  we  arrive  at  the  equation 


(3.22) 


-  2(7o  -  1) 


-«2 


>1-  72 


a^“  7i(7i-  1)(72  -  l) 

and  for  the  roots  of  equations  (3»i9)  and  (3«20)  we  obtain  at  the 
same  tine  the  expressions 


23) 


(3.24) 


It  is  obvious  that  tb*  corresponding  solution  determines  the  point 
located  on  the  parabola 

Zg  -  (Vg  -  if.  (3.25) 

If  the  inequality  (3*22)  exists,  then  two  solutions  are  obtained  for 

which  the  corresponding  points  M  and  N  (Fig.  3)  are  located  on  different 

sides  of  the  parabola  (3*25). 

It  is  not  difficult  to  satisfy  oneself  that  if  the  ratio 
2 

(qi  -  assigned  in  such  a  way  that  the  Inequality  (3.l8)  is 

p  p 

satisfied,  then  for  different  values  of  c  /a^  the  corresponding  points 
M  and  N  will  be  located  on  a  certain  third  order  curve,  while  eqtiation 
(3.23)  is  satisfied  in  the  points  of  the  parabola  (3*25),  and  the  quantity 
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2  2 

c  /a^  attains  an  extremum.  The  largest  root  of  equation  (3.25) 
corresponds  to  a  minimum,  and  the  smallest  to  a  maximian  of  the  ratio 


If  in  relation  (3.21)  there  is  the  sign  of  equality,  then  equation 
(3.25)  has  the  form 


(3.26) 


It  is  obvious  t^at  in  this  case  the  velocity  c^  is  equal  to  the 
velocity  of  sound  in  a  gas  vith  the  adiabatic  coefficient  7^.  From 
formulas  (3.24)  for  the  corresponding  values  of  and  we  find 
Vg  =>  0,  and  Zg  1. 


If  in  the  relation  (3.21)  there  is  a  sign  of  inequal ity,  then  the 
equation?  (5.24)  have  two  roots  c^  and  Cg  that  satisfy  the  inequality 


All  the  reasonings  of  this  paragraph  are  applicable  to  motions  with 


plane  waves  and  to  motion  with  cylindrical  and  spherical  synanetry. 
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4.  Let  us  now  consider  the  q,ue«tlon  of  determining  ell  the 
Bolutione  when  m  -  2,  n  -  -2,  s  •  0,  k  -5,  eisd  V  »  5.  Let  us 

assx&oc 


a  -  Pt  ,  b 


.  o  2  t 

^  “  *1  "2  " 


where  and  are  certain  selected  valvuss  of  pressure  and  density. 
Equations  (2.5)  and  (2.6)  in  the  case  under  consideration  assvime 


the  form 


d*  _  „  Z  *-(v-i)(rv-  1) 


1  d  In  \ 

7  “(iV 


5*-(V-  1)^ 

vr5z-(v- 1)' 


i^.l) 


(^.2) 


Moreover j  when  V  and  the  last  of  equations  (2.5)  is  Integrated, 
we  obtain  the  folloifing  connection  between  P,  R,  and 


.7-1  * 


(4.5) 


where  ^  is  a  constant  of  integration. 

Using  the  relation  (4.5)  and  the  designation  yP/R  -  t,  we  can  write 


fomiolas  (1.1)  in  the  form 


pop 


(4.4) 


In  regions  of  continuous  lootion  the  quantity  p/p^  of  the  various 
fluid  particles  will  be  identical,  since 
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JL. 


Pi 


1 


(1^.5) 


It  is  obvious  that  the  abstract  constant  p  is  closely  connected  with 
the  value  of  entropy.  If,  at  a  certain  point  in  the  region  of  continuous 
notion,  pressiu'e  and  density  have  the  values  and  p^,  then  it  is 
obvious  that  in  this  region  p  »  1. 

The  solution  eiqpressed  by  formulas  (4.4-)  depends  on  four  arbitrary 
constants.  Two  of  them  are  and  p^P;  two  others  vlU  enter  during 
the  integration  of  the  differential  equations  (4.1)  and  (4.2). 

The  main  difficulty  lies  in  the  integration  of  equation  (4.1). 

After  the  integration  of  this  equation,  the  parameter  X.  is  determined 
as  a  function  of  V  from  equation  (4.2)  by  means  of  simple  quadrature. 

Let  us  consider  the  behavior  of  the  Integral  curves  of  equation 
(4.1)  in  the  aV  plane.  This  equation  has  five  special  points: 

0(0,0),  A(1,0),  B(t*,7*),  C(0,1),  D(0D,  +a0), 

while 


»* 


(57- 1) 


5  ^ 


The  points  0,  A,  and  C  are  nodes  and  the  points  B  and  0  are  saddles. 
The  general  form  of  the  integral  oturves  in  the  half-plauie  z  >  0  is  shawn 
in  Fig.  4. 


Fig.  4 
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In  the  points  of  the  parabola  c  <■  (V-  l)(7V-  l),  the  Integral 
curves  have  tangents  parallel  to  the  V-axls;  in  the  points  of  the 
parabola  (Y~  >  they  have  tcuigents  parallel  to  the  t-axis. 

Vear  the  point  0«  the  solutions  of  equatiens  (^.l)  and  (^.2) 
have  the  form 

.2 


g  2/2 

where  >  0  and  kg  ^  0  sure  aihitrary  constants.  \  ■  a^^  t  /r  -  0 
corresponds  to  the  point  0;  conscQ.uently,  the  infinitely  irenote  point 
r  •  OD  corresponds  to  the  point  0  in  the  space  of  the  notion  of  the 
fluid. 

At  infinity  for  velocity,  density,  and  pressure  the  values  v^, 

Pq  and  p^  are  detemined  by  the  fozmulaa 


F* 


If  p 


(y-i) 


o  Pi'  Po  ■  Pi'  P  -  1,  kj^  -  kg 


In  the  B,  V  plane  near  the  point  0,  the  Integral  curves  have  the 
form  of  parabolas.  As  the  velocity  increases  to  Infinity  the  corres¬ 
ponding  parabolas  approach  the  V  suds, 
upward  along  the  right-hand  bremch  of  the  parabola  corresponds  to  the 


If  v^  >  0,  then  the  shifting 


approach  from  infinity  to  the  center  of  synmetry;  if  v^  <  0,  then  the 
shifting  vqpvard  along  the  left-hand  branch  of  the  parabola  corresponds 
to  the  approach  from  infinity  to  the  center  of  syonetry. 

Xquatlons  (k.l)  and  (^.2)  have  an  exact  solution: 

^2 


*  -  0, 


V  - 
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Th«  condition*  of  motion  In  vhlch  either  p  »  0  or  p  •  co  corree- 
pond  to  this  solution. 

Moreover,  eqlustions  (4.1)  and  (4,2)  have  the  solution 

y,  -  0,  a  »  kj\  .  (4.9) 

It  is  obvious  that  this  solution  corresponds  to  the  state  of  rest 

V  -  0,  p  »  p  -  Pj^pkj^  (4.10) 

If  p^  and  are  the  density  and  press\ire  correspoxuling  to  the 

state  of  rest,  then  0  «  1,  k^  •>  1. 

I^e  integrated  straight  line  satisfies  the  state  of  rest.  Ihe 

2  2/2 

moving  points  of  the  fluid  for  which  k  «  a^  t  /r  »  constant  corresitoztd 
to  the  fixed  p9ints  of  this  strai^t  line. 

It  is  not  difficult  to  see  that  to  the  point  z  «  1,  V  ■  0  there 
corresponds  a  state  of  motion  which  propagates  in  the  fluid  with  a 
velocity  eq^ual  to  the  velocity  of  sound  satisfying  this  state. 

Indeed,  from  fomulas  (4,9)  end  (4.10)  when  z  ■  1  we  obtain 


When  z  <  1,  the  velocity  of  propagation  of  the  corresponding  states 
in  the  quiescent  fluid  is  greater  than  the  velocity  of  sound,  and  lAien 
z  >  1  it  is  less  than  the  velocity  of  sound. 

Let  us  consider  the  behavior  of  the  integral  curves  near  the  point 
A  for  which  z  ■  1,  V  «  0.  Hear  this  point  the  equation  of  the  first 
appixoxlnation  for  equation  (4.1)  has  the  form 


dz  _  z-1  +  (1+  7)V 

5^ - H 


(4.11) 
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The  general  Integral  of  equation  (4.11)  is  represented  by  the 


formula 

z-1  -  (TV  +  (1+  7)V  In  V,  <4.12) 

vhere  C  Is  a  constant  of  integration. 

It  Is  obvious  that  at  the  point  A  there  is  a  nodsj  and  that  the 
Integral  cvirves  converge  to  the  point  A,  tangent  to  the  z  axis. 

It  is  not  difficult  to  see  that  in  the  motion  along  the  integral 
curve,  the  point  A  is  reached  for  a  finite  value  of  ths  parameter 
since  at  the  point  A  ve  have 


During  the  motion  to  the  point  B  along  the  integral  curves  corres¬ 
ponding  to  the  plus  sign,  the  parameter  \  tends  toward  during  motion 

to  the  point  B  along  the  integral  curves  corresponding  to  the  minus  sign. 


the  parameter  X  tends  tovnzrd  zero. 


A  special  solution  satisfies  the  point  B: 


(4.15) 
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The  solution  (4.13)  is  b.  pairtlcular  case  of  the  d«p«nilence  of 
the  general  solution  on  the  one-dliaenalonal  constant  vhlch  vc  examined 
In  paragraphs  1  and  2.  Formulas  (l.l)  and  (2.4)  become  fozmulauB 
(4,15)  when 


k  «  ,  8  ■  .  (4<,l6) 

In  the  solution  of  (4.15),  when  t  «  0,  r^O,  we  have  v  -  aO  , 

p  ■*  at),  and  p  «  os ;  when  t  >  0,  the  infinite  values  are  maintained  only 

when  r  ■  op;  when  r  =■  0,  we  have  v  •  0,  p  •  0,  p  ■  0. 

Let  ua  consider,  finally,  the  special  point  C.  Hear  the  point 

V  «  1,  z  «  0  for  sufficiently  small  ^  ^  (V-1)  ,  and  the  continuation 

of  the  integral  curves  to  the  right  auid  left  lead  to  the  Intersection 

with  the  parabola  *  "  y  (V-l)  .  In  the  points  of  this  parabola  we  have 

(dz/ov)  -00;  during  fuller  continuation  the  integral  curves  tend  toward 

1  2 

the  point  C,  remaining  below  the  parabola  t  -  ^V-1)  .  Consequently 
the  point  C  is  a  m^e. 

It  is  not  difficult  to  show  that  for  all  values  of  7  >  1  the 
following  Inequality  exists: 


lim  - 5 —  <  i  . 

V-*!  (V-1)  ^ 

From  this  it  foU-owe  that  the  derivative  d  In  \/dV  has  a 
finite  value  when  V  “  1.  Consequently,  during  the  motion  along 
the  integral  curve  we  arrive  at  the  point  C  with  a  finite  value 
of  the  parameter  The  point  C  may  correspond  to  a  stationary 
discontinuity  or  to  a  boundary^  which  is  shifting  together  with 
the  fluid  particles. 


(l^.l?) 
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At  the  point  C  we  have  2  •»  yf/'B.  -  0,  and  from  this  and  fran 
(^.3)  it  follows  that  for  all  solutions  different  from  t  »  0,  dur¬ 
ing  the  q[>proach  to  the  point  C  pressure  and  density  tend  toward 
zero. 

Motion  along  the  Integral  curve  on  the  side  of  increase  in 
parameter  X.  corz^esponds  to  a  decrease  in  r.  When  0  <  r  <  03  ,  we 
have  aO  <  \  <  0. 

By  meauiB  of  equation  (4.2)  it  is  easy  to  show  the  directions 
of  increase  in  parameter  X  along  the  Integral  carves.  In  Fig.  4 
these  directions  are  shown  by  arrows. 

On  the  parabola  z  •  (V-l)  ,  the  parameter  X  reaches  a  raaximvm 
when  V  <  0  and  V  >  1  and  a  minimal  when  0  <  V  <  1. 

It  is  obvious  that  the  continvious  transition  along  the  integral, 
curve  through  the  points  of  this  parabola,  exclxxling  the  points  A 
and  C,  is  physicsLlly  inadzaissible,  since  it  leads  to  non-ainglevaluedness. 

In  the  center  of  symmetry  of  the  stream,  r  »  0  and  X  «  cD  . 

In  the  zV  plane,  the  motion  of  the  fluid  continues  to  the  center 
r  ■  0  along  the  integral  straight  line  V  -  0,  z-^+oO,  satisfying 
the  state  of  rest;  along  the  atrai^t  line  z  -  0,  V->+ oO,  for  which 
either  the  density  is  infinite  or  the  pressure  is  equal  to  zero, 
for  the  special  solution  (4.15) #  i.e.,  for  the  point  B,  and  finally, 
for  the  motion  which  satlaflea  the  Integral  curves  OB  and  DB,  resting 
on  a  special  solution  for  the  point  B. 
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At  «n  infinitely  remote  point  of  the  etreani  we  have  r  fltt  , 

\  ■  0.  In  the  plane  only  the  special  points  0  and  B  can  corres¬ 
pond  to  the  infinitely  remote  point  of  the  stream.  For  the  point  B, 
the  velocity  is  infinite  at  Infinity.  Pbr  solutions  starting  from 
the  point  0,  the  velocity  of  the  fluid  at  infinity  can  have  any 
value. 

2*  Based  on  the  general  analysis,  given  above,  of  the  conditions 
at  the  jung>s  and  the  behavior  of  the  integral  curves  of  equations 
(4.1)  and  (4„2)  in  the  plauie,  it  is  not  difficidt  to  construct 
solutions  of  the  problems  pointed  out  in  Section  1  (page  2 ) . 

Problem  1°.  The  form  of  the  integral  curve  when  v^  <0  is 
shown  in  Fig.  5.  The  initial  data  for'Sae  integral  curve  are  determin¬ 
ed  by  the  values  •  v^,  ■  Pq»  •nd  p^  •  by  means  of 

'formulae  (4.7)  and  (4.8). 

At  the  point  of  the  intersection  of  the  parabola 
X  •  (l-V)  fl+(7-l)V/^wlth  the  integral  curve,  a  strong  Jump  aurises, 
beyond  which  the  fluid  is  at  rest  (point  Hg) .  The  velocity  of  pro¬ 
pagation  of  the  Jx*i?)  Of  which  is  equed  to  the  velocity  ox'  growth  of 
the  radius  of  the  spherical  nucleus  with  the  stationary  gas,  is 
determined  by  the  value  of  the  parameter  \  -  XPj^/(p^c  ),  satisfying 
the  points  and  Ng. 

If  Vj^  >0,  then  the  corresponding  integral  curves  are  represented 
in  fig.  6.  For  sufficiently  smalT  vadues  of  v^,  integral  curves  of 
the  type  OA  ere  obtained,  the  motion  being  continuous  for  edl  values 
of  r.  Fear  the  center  of  syta&etry  a  q;>berical  nucleus  of  quiescent  gas 
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forms  that  trsmsforms  by  a  weak  diecontinuity  into  motion,  and  the 

point  A  corresponds  to  the  weak  discontinuity.  At  a  certain  initial 

« 

velocity  ^  integral  curve  OB  is  formed  in  th€  plane. 

The  point  B  corresponds  to  the  center  of  syjanetry;  in  this  case  the 

velocity  of  the  gas  will  be  equal  to  zero  only  at  one  point  -  the 

* 

center  of  symmetry.  If  "'^3^  >  then  the  integral  curves  converge 

toward  the  point  C,  at  which  the  parameter  \  has  a  certain  finite, 

* 

nonzero,  constant  value  k  ,  and  pressure  and  density  become  zero. 

* 

Consequently,  at  the  initial  velocities  v^  >  v^  in  the  gas  a  vacuum 
is  formed,  expanding  with  constant  velocity  determined  by  the  value 
\  . 

Problem  2^.  Let  us  designate  by  U  the  velocity  of  a  spherical 
piston. which  is  equal  to  the  velocity  of  the  gas  particles  on  the 

sphere  the  pressure  inside  this  sphere  is  constant  and  equal  to 

* 

P  • 

We  designate  hy  and  the  density  and  pressure  in  the 
quiescent  fluid.  In  the  zV  plane  a  certain  point  of  the  straight 
line  V  =  1  corresponds  to  the  bounding  sphere  3^.  Integral  curves, 
starting  frcan  points  of  the  straight  line  V  ■  1,  intersect  the 
pai'abola  z  -  (V-l)  .  Therefore,  continuation  of  motion  to  the  point 
0,  which  cori*e8ponds  to  an  infinitely  remote  point,  is  possible  only 
by  a  Junq). 


Fig.  7 


Assuming  that  the  motion  of  the  fluid  starts  from  the  state 
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of  rest,  it  follows  that  on  the  external  side  of  the  «  0, 

z^^O,  i.e. ,  the  corresponding  point  is  located  on  the  z-axis,  whence 
it  follows  that  the  point  Kg(zg,  is  detemnined  hy  the  intersection 

of  the  integral  curve  (|*ig.  7)  with  the  parabola 


z  -  (1-V)(1  +  V). 

« 

On  the  straight  line  V  <•  1,  the  point  M 

sphere  S  and  the  value 
o 


corresponds  to  the 


* 

\ 


>  K 


7Pl 


so  that  c  >  U. 

2  2  *  * 

When  c  — ^  ,  the  values  z  and  K  tend  to  infinity,  and  the 

velocity  U  tends  to  zeix3. 

As  a  result  of  numerical  calculations,  we  have  determined  the 
dependence  of  the  ratios  U/c,  p  /p^,  p*/p^  on  the  quantity  a^/c. 

In  Fig.  2  these  dependencies  are  represented  in  the  form  of  dashed 
curves;  coinparison  of  the  dashed  curves  with  the  solid  carves  gives 
an  increase  in  the  ratios  under  consideration  at  the  transition 
from  the  jvmp  to  the  piston  (sphere  S^)* 

The  examined  motion  may  be  continued  inside  the  sphere  S^. 

During  the  continuation  the  5)here  may  become  a  surface  of 
stationary  discontinuity  (only  density  is  discontinuous;  velocity 
and  pressure  are  continuous).  On  the  straight  line  V  =*  1,  the  integral 
curve  z  =■  z(V)  may  have  a  point  of  discontinuity.  The  magnitude  of 
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#  I 

discontinuity  (location  of  the  points  M  and  M  )  is  determined  by 
the  ratio  of  the  denaitiee  on  the  different  sides  of  the  sphere 
S^,  which  may  he  given  nrhiti^rlly. 

Prom  the  condition  of  ainglevaliJiednesB,  it  follows  that  from 

I 

the  point  M  it  is  possible  to  continue  the  integral  curve  to  the 

riglit  (see  Fig.  7)  only  to  the  point  of  the  Intersection  with 

the  parabola  z  (V-l)  .  After  the  intersection  of  the  integral 

2 

Cuzv*  with  the  parabola  z  =  27(7-1)  /(7-I)  at  the  point  E,  further 

continuation  of  motion  may  be  realized  by  a  Jump  from  a  certain 

point  of  the  segment  ED  to  the  corresponding  point  Hj^.  From  the 

point  the  integral  curve  may  be  continued  only  when  the  conditions 

of  singlevaluodness  to  the  intersection  of  the  parabola  z  =>  (V-l) 

at  the  point  Dg  are  satisfied.  According  to  the  theorem  of  Tsemplen 

1 

no  other  Juiups  can  occur.  The  velocities  of  propagation  c,  U,  o  , 

<  *  • 

and  U  correspond  to  the  points  5^;  MM  ;  and  Dg.  The 

I  I 

oocrectness  of  the  inequalities  e>U>c  >U  is  obvious. 

Consequently,  the  motion  con  be  continued  inside  S^,  but  this 
continuation  can  be  realized  only  up  to  a  certain  sphere  S,  the 

f 

radius  of  which  increases  according  to  the  law  r  *  U  t.  On  the 
surface  of  the  sphere  S,  the  density,  pressizre,  and  velocities  of 
the  fluid  (greater  than  U  )  are  constant. 

Problem  On  the  front  of  the  spherical  detonation  we  have 

2/2  2/2 
\  w  a^^  /e  «  constant.  The  point  ^2.  ^  ,  V  ■  0  corresponds  to 

the  external  side  of  the  Jub^.  The  value  of  all  the  characteristics 

on  the  internal  side  of  the  jun^)  will  be  determined  when  the  follow- 
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ing  are  assigned:  the  density  and  pressure  in  the  quiescent 
mixture,  the  heat  of  reaction  “  ^2  ^^at  satisfies  the  ineqtiality 
(5.21),  and  the  velocity  of  propagation  of  the  shock  vave  c.  Prom 
the  conditions  at  the  two  solutions  corresponding  to  the  points 

M  and  N  (see  Hg.  3)  are  obtained. 

Continuation  of  notion  from  the  point  M  Is  accon^jlished  by  the 
integral  curve  MA  (see  Fig.  8).  The  point  A  corresponds  to  the 
boundary  of  the  expanding  nuc2.eu8  in  which  the  gas  -  the  product  of 
the  chemical  reaction  acconq^suaying  the  detonation  -  is  at  rest. 


Continuation  of  motion  from  the  point  H  is  accomplished  by 
the  integral  curve  NM  .  It  is  obvious  that  it  is  Impossible  to 
continue  this  motion  to  the  center  of  synmetry.  Motions  of  such 
character  may  be  excited  by  metuas  of  the  additionsJ.  spherical 
piston  which  we  considered  in  Problem  2°. 

The  velocity  of  propagation  of  the  detonation  at  a  given  value 
q^  -  qg  is  a  minimum  when  the  points  M  and  N  converge  into  a  single 

p 

point  Q  locatdd  on  the  parabola  z  -  (V-l)  .  In  this  critical  case 
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the  motion  can  be  continued  either  to  the  right  at  the  point 
or  to  the  left  to  the  point  A.  In  the  first  case,  pressures  greater 
than  at  the  Jud^  are  obtained  In  the  Internal  region  and  on  the 
piston;  in  the  second  case,  beyond  the  Jun?)  a  rarefaction  vave  ^ 
arises  in  vhlch  the  pressures  diminish  from  the  point  Q  vip  to  the 
point  A. 

Various  autlaore  paid  much  attention  to  the  problem  of  detennln- 
Ing  the  detonation  velocity.  We  con  point  out  that  all  hydrodynamic 
conditions  are  satisfied  if  a  Jump  takes  place  at  the  point  M  and 
the  motion  continues  along  the  integral  curve  MA. 

Experimental  data  show  that  there  is  a  regime  of  detonation 
coi’responding  to  a  minimum  velocity  satisfying  the  point  Q  (see  Fig. 

8)j  with  subsequent  change  of  chamcteristics  along  the  integral 
curve  QA.  substantiation  of  this  fact,  Zeldovlch  drew  vpon  scane 

reasonings  connected  vkth  the  chemical  mechanism  of  the  reaction  in 
the  shock-vave  front. 

It  is  obvious  that  at  any  moment  each  motion  of  the  type  under 
consideration  is  mechanically  similar  to  ajiy  state,  and,  in  particvilar, 
to  the  state  for  the  moment  of  time,  as  close  as  desired  to  the 
initial  moment  of  detonation  at  the  center  of  symmetry.  This  circum¬ 
stance  con  serve  as  a  reason  for  the  choice  of  velocity  of  detonation 
from  the  data  on  the  method  of  excitation  in  the  center  of  symmetry. 

If  the  detonation  is  initiated  in  the  center  by  external  pressure, 

^For  example,  see  Zeldovlch  f^l  and  Qrlb  M  • 
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then  a  point  of  type  H  rauut  correspond  to  the  Junqp.  On  the  other 
hand,  the  demand  that  the  velocity  of  the  products  of  reaction  at 
the  center  of  synnetry  should  equal  zero  leads  to  solutions  in 
•which  a  point  of  type  M  corresponds  to  a  Jump.  Simultaneous  ful¬ 
filling  of  these  "two  conditions  ceui  he  satisfied  only  for  points 
of  type  Q  and  an  Integral  curve  of  the  type  These  considerations 

may  be  considered  as  a  certain  basis  for  the  fact  that  in  practice 
a  regime  of  jumps  corresponding  to  points  of  the  type  Q  m\ist  be 
obtained. 

Received  by  the  editor 
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